Spin dependent structure function gi at 

low X and low 



B. Badelek J. KiRYLUK ^ and J. Kwiecinski 



Department of Physics, Uppsala University, P.O.Box 530, 751 21 Uppsala, Sweden 
^ Institute of Experimental Physics, Warsaw University, Hoza 69, 00-681 Warsaw, 
Poland 

Department of Theoretical Physics, H. Niewodniczanski Institute of Nuclear Physics, 
Radzikowskiego 152, 31-342 Cracow, Poland 

e-mail addresses: badelek@fuw.edu.pl, kiryluk@fuw.edu.pl, jkwiecin@solaris.ifj.edu.pl 

Abstract 

Theoretical description of the spin dependent structure function gi{x,Q'^) 
in the region of low values of x and is presented. It contains the Vector 
Meson Dominance contribution and the QCD improved parton model suitably 
extended to the low domain. Theoretical predictions are compared with the 
recent experimental data in the low x, low region. 

1. Introduction 

Measurements of polarised deep inelastic lepton-nucleon scattering have determined 
the cross section asymmetries Ai and spin dependent structure functions gi of the 
proton, deuteron and neutron in a wide kinematic range of and x.Q This allowed a 
verification of sum rules, like e.g. the Bjorken sum rule which is a fundamental rela- 
tion in the QCD, and the Ellis- Jaffe sum rules. Evaluation of the sum rules requires 
knowledge of the structure functions gi over the entire region of x as well as their evo- 
lution to a common value of Q^. Since the experimentally accessible x range is limited, 
extrapolations to a; = and x = 1 are necessary. Of these the former is critical since 
the small x behaviour of gi{x) is theoretically not well established and the relevant 
contribution to the sum rules' integral may in principle be large. 



Theoretical predictions for the structure function gi over a full range of x are even 

more interesting than for its first moment, especially at low x, i.e. at high parton 

densities, where the new dynamical mechanisms may be revealed. Theoretical and 

experimental studies at low x in the polarised case are thus awaited for. A possible 

""^Here, as usual, x = /{2pq) where = —q^ with q and p denoting the four momentum transfer 
between leptons and the four momentum of the nucleon respectively. 
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future polarising the proton beam at HERA would be a milestone in this field. 

In the fixed target experiments the low values of x are reached by lowering at the 
same time the values of Q^. Theoretical analysis of these data therefore requires a suit- 
able extrapolation of the structure function to the low region. Low phenomena 
and in particular a transition from the perturbative (large Q^) to the nonperturbative 
(low Q^, including (5^=0) region is actively investigated in the spin insensitive exper- 
iments. In spite of a wealth of data and of a wide spectrum of ideas this field is still 
a major challenge in high energy physics |]T|. Among the spin sensitive experiments 
the only available low data are from the E143 experiment at SLAC (moderate 
X and low Q'^) and now also from the SMC at CERN [§, ^ (low x and low Q^). In 
the low region one can expect that dynamical mechanisms, like the Vector Meson 
Dominance (VMD), can play an important role. For large the VMD contribution to 
gi gives a power correction term and can usually be neglected. Moreover, the partonic 
contribution to gi which controls the structure functions in the deep inelastic domain 
has to be suitably extended in order to get its extrapolation to the low region. 
The latter component will be expressed in terms of the unintegrated (spin dependent) 
parton distributions and we show that the corresponding representation of gi can be 
easily extrapolated to the low region. The main purpose of our paper is therefore 
to construct the structure function gi^x.Q"^) which would include the VMD and the 
(QCD improved) parton model contributions. 

The content of the paper is as follows: in the next Section we present the data on 
gi and comment on the Regge model predictions for gi which are often being used for 
x = extrapolations. In Sec. 3 we briefiy present a formalism describing gi in terms 
of the unintegrated spin dependent parton distributions, incorporating the leading or- 
der Altarelli-Parisi evolution and the double logarithmic ln^(l/x) resummation at low 
X. In Sec. 4 we discuss the Vector Meson Dominance part of the gi which has to be 
included in that region since, as it has already been pointed out above, for the fixed 
target experiments low values of x are correlated with the low values of Q^. Numeri- 
cal results are also presented there. Finally, in Sec. 5 we give a summary of our analysis. 



2. The gi data 

Several experiments contributed to the spin structure function gi measurements on 
different targets and over different kinematic intervals. As a result, for proton and 
deuteron, gi was measured for 0.00006 < a; < 0.8 by the EMC |], SMC |, |, E143 0, 
E155 [^] and HERMES [^. For neutron, gi was measured for 0.014 < a; < 0.8 by the 
E142 H, E154 H and HERMES [|lOl. A summary of xgl''^{x) data at the measured 
values is presented in Fig.l. 
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For the SMC data, (x) = 0.0001 corresponds to (Q^) = o.02 GeV^. In other expe- 
riments Qi was measured with high statistical accuracy for x > 0.01 and for > 1 
GeV^ only.0 We do not present as there are no direct measurements for x < 0.01 
i.e. in the low x region. 



The lowest x and region was explored by the SMC due to a high energy 
of the muon beam and implementation of a dedicated low x trigger. The results 
of the SMC presented in Fig.l come from two different analyses [Q, ^ which join 
at X ~0.002. It should be noted that a direct result of the measurements is the 
virtual photon-nucleon asymmetry, Ai. To get the gi one has to use the relation 
gi = Ai- Fi= Ai- F2/[2x{l + R)], where F2 = Fl + Ft, R = Fl/Ft and Ft = 2xFi 
with Fl and Ft denoting the unpolarised nucleon structure functions corresponding 
to longitudinal and transverse polarisations of the virtual photon respectively. Unfor- 
tunately there have been no direct measurements of F2 and R in the kinematic region 
of the low X and low gMC data, i.e. for 0.00006 < x < 0.003 and 0.01 < < 1 
GeV^. Thus the SMC used the model for the F2 and a parametrisation of Ref.|]12 
for R so their results for gi are model-dependent. 



The new low x data of the SMC include the kinematic region where W"^ = {p + q)"^ 
is high, W'^ > 100 GeV^ and much larger than Q^. Thus one should expect that the 
Regge model should be applicable there. According to the Regge model, gi{x,Q^) ~ 
^-a for X ^ and fixed Q^, where a denotes the intercept of the Regge pole trajectory 
corresponding to axial vector mesons. It is expected that a ~ for both / = and 

1 = 1 trajectories, |13|. This behaviour of gi should go smoothly to the W'^°' depen- 
dence for Q'^ — i> 0. Other considerations related to the Regge theory predict gi ~lnx, 
[ll4| , while the model based on exchange of two nonperturbative gluons gives gi ~ 

2 ln(l/a;)-l, |T5[. A perverse behaviour, gi ~l/(xln^x), recalled in [|T^, is not valid for 

91, m- 



In the kinematic range of the SMC data chan ges very little: from about 100 
GeV^ at X = 0.1 to about 220 GeV^ at x = 0.0001, contrary to a quite strong change of 

(from about 20 GeV^ to about 0.01 GeV^ respectively). This means that the new 
SMC measurements cannot test the Regge behaviour of gi through the x dependence 
of the latter, without additional assumptions about the Q"^ dependence of gi. A model 
which allows extrapolation of gi to the low region is described in the next Section. 

^The E143 measured the asymmetry Af^'" for 0.024 <x < 0.205 and 0.31 < < 1 GeV^ but gi 
was not extracted from those data. 
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3. Partonic contribution to gi 



In the region of large values of the spin depedent structure functions are described 
by the QCD improved parton model 0]. In this model gi = gf""^^, where gf""^^ is 
related in a standard way to the polarised quark and antiquark distributions Aqi and 
Aqi corresponding to the quark (antiquark) flavour i: 

gr\x,Q') = l E e^[Aq^{x,Q') + Aq,{x,Q')]. (1) 

i=u,d,s 

In what follows we assume Aq^ = Aq^ and set the number of flavours equal 3. 

In perturbative QCD the structure function g^'^^'^ is controlled at low x by the double 
logarithmic ln^(l/x) contributions i.e. by those terms of the perturbative expansion 
which correspond to the powers of ln^(l/x) at each order of the expansion It is 



convenient to discuss the ln^(l/x) resummation using the formalism of the unintegrated 
(spin dependent) parton distributions fj{x',k'^) (j = Uv,dy,u,d,s,g) where /c^ is the 
transverse momentum squared of the parton j and x' the longitudinal momentum 
fraction of the parent nucleon carried by a parton []T^ |2T|. The conventional 
(integrated) distributions Apj{x,Q^) (i.e. Ag„ = Ap„^ + Ap^, Ag„ = Ap^ etc. for 
quarks, antiquarks and gluons) are related in the following way to the unintegrated 
distributions fj{x',k'^): 

Ap,{x,Q') = Ap'^{x)+J^^^ -^a^' = <^+gl)^k') (2) 

Here Ap^-[x) denote the nonperturbative parts of the of the distributions, corresponding 
to fc^ < /cq the parameter is the infrared cut-off {k^ ~1 GeV^). They are treated 
semiphenomenologically and parametrised in the form used in Refs |jT^, ^ : 

Ap]{x) = C,{l-xf^ (3) 



In Eq.(^ we assumed ?7^l„ = rjd^ =3, r]u = r]s = 7 and rjg = 5. We also used kQ =1 
GeV^. The normalisation constants Cj were determined by imposing the Bjorken sum 
rule for — A(i° and by requiring that the flrst moments of all other distributions 
are the same as those determined from the QCD analysis of [E^ . 



The unintegrated distributions fj{x', k"^) are the solutions of the integral equations 
[]T^ |2T| which embody both the LO Altarelli-Parisi evolution and the double 



ln^(l/x') resummation at small x' . These equations combined with equations (P and 
@) lead to approximate x~'^ behaviour of the g^""^* in the x — > limit, with A ~ 0.3 
and A ~ 1 for the nonsinglet and singlet parts respectively which is more singular at 
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low X than that generated by the (nonperturbative) Regge pole exchanges [|. The dou- 
ble ln^(l/x) effects are presumably not important in the W"^ range of the fixed target 
experiments (cf. Fig. 2 in and Fig. 6 in [^) but they significantly affect gi in the 
low X region which may be probed at the polarised HERA, [|l^, 20, 



Va. However the 



formalism based on the unintegrated distributions employed here is very suitable for 
extrapolating gi to the region of low at fixed W'^ [|19| . 



Formulae (|l]) and @ define partonic contribution to the structure function gi{x^ Q"^). 
Since x{l + k'^/Q'^) k'^/W'^ for — * in the integrand in Eq. (j^) and since k'^ > k^ 
there, the gi"'^^{x, Q"^) defined by Eqs (|l|) and (0) can be smoothly extrapolated to the 
low region, including = 0. In that limit, the gi should be a finite function of 
W^, free from any kinematical singularities or zeros. The extrapolation, valid for fixed 
and large W'^, can thus be done for the gi"'^^{x, Q"^) given by Eqs (JlD and (||) provided 
that nonperturbative parts of the parton distributions Apj(x) are free from kinematical 
singularities at x = 0, as in the parametrisations defined by Eq. (H). If Ap°(x) contain 
kinematical singularities at x = then one may replace Ap°(x) with Ap°(x) where 
X = x{l + kl/Q'^) and leave remaining parts of the calculation unchanged. After this 
simple rearrangement the structure function gi°'^*{x, Q"^) can be extrapolated to the 
low region (for fixed W"^) including the point = 0. Possibility of extrapolation 
to = is an important property of the formalism based on the unintegrated parton 
distributions. 



We solved equations for the functions fi{x\k'^) [T^, |20|, ^ and calculated the 
gi"'^^{x,Q'^) from Eqs (|5) and (|^) using the parametrisation (^. To be precise we 
solved equations which resummed only the ladder diagrams contributions in that part 
which corresponded to the double ln^(l/x) resummation but this approximmation was 
completely adequate for the values of which are relevant for the fixed target ex- 
periments. Let us also remind that equations for the functions |]19|, ^ 
combined with equations (|l],|D are a generalisation of the LO QCD evolution equations 
p3| for polarised parton densities and for moderately small and large values of x are 
equivalent to these equations. 



As a consequence gi"'^^ calculated at x and values of the SMC measurement 
gives a reasonable description of the SMC data on gi''^{x,Q'^), cf. Fig.l (it does not 
reproduce at the same time other measurements equally well due to differences in 
values between the experiments). For the sake of the comparison the calculated gi""^^ 
was extrapolated to low values of since all the data with x < 0.001 have < 1 

To be precise the singular x~'^ behaviour with A ^ 1 for singlet and gluon spin dependent 
distributions does hold in the approximation when only the ladder diagrams are retained pO{ . Com- 
plete double logarithmic ln?{l/x) resummation which includes also the non-ladder bremmstrahlung 
diagrams generates less singular behaviour of these distributions j2l). 
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GeV^. However the (extrapolated) g^""^*^ may not be the only contribution to gi in the 
low domain. 

4. Vector Meson Dominance contribution to gi 

One expects that in the low region an important role may be played by the VMD 
mechanism. The structure function should thus be represented by the sum of the 
partonic and VMD contributions, i.e. 

9i{x, Q') = gX''^{x, Q') + gr\x, Q') (4) 

The VMD contribution to gi{x, Q"^) can be written as: 

In this formula the constants 7^ are determined from the leptonic widths of the vec- 



tor mesons [24| and denotes the mass of the vector meson v. The cross sections 



AaviW^) are for high energy given as the following combinations of the spin de- 
pendent total cross sections: 

Cri/2 - (73/2 , s 

^(^v = -^-^ (6) 

where o"i/2 and cr3/2 correspond to the total vector meson - nucleon cross sections 
with the projections of the total spin on the vector meson momentum equal 1/2 and 
3/2 respectively [^. Unfortunately the cross-sections Aa^ are unknown. In order to 
estimate the VMD contribution, gY'^^{x,Q'^), we assume that the cross sections Acr^ 
are proportional to the appropriate combinations of the nonperturbative contributions 
Ap°(x), defined by Eq.(|^), to the polarised quark and antiquark distributions. For the 
proton we assume: 



E 



pq 

47r v^,co IviQ^ + 



C 



^ [Aulix) + 2Au\x)) + \ [Adlix) + 2Au\x\ 



4 

ml 



(g2 + m 



2\2 



(7) 



pq 



mtAasv „2 . ml 



47r7|(Q2 + ^2)2 9 W(g2 



Ct^^Ax) ,^, ' (8) 



m 



where AvP{x) = Ap^{x) etc. All distributions are parton distributions in the proton. 
The distributions Ap°(x), Eq. @, behave as x° for x — > 0. As a result the cross 
sections AcTy behave as l/W"^ at large that corresponds to the assumption that the 
corresponding Regge trajectories have their intercepts equal to zero. We include ex- 
act X dependence of the nonperturbative (spin dependent) parton distributions Ap°(x) 
and not only their (constant) x Q limits, Cj. This gives an extension of the VMD 
model to the region of moderately small values of x. Formally this means that we allow 
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additional dependence of the cross-sections Acr^ in terms which are non-leading in 
the large W"^ limit, i.e. vanish faster than l/W"^. 

We shall vary the parameter C in Eqs and (^) and analyse a dependence of the 
structure function gi{x, Q^) upon the value of this parameter. It should be noted that 
the VMD part of gi vanishes at large as l/Q^ (contrary to the g^'^^^ which scales 
modulo logarithmic corrections) but it may be a dominant contribution at (very) low 

as it is the case for the unpolarised structure functions. For low we expect a 
dominance of the VMD part of gi. In analogy with the unpolarised case we expect 
that it should exhaust about 80 % of the total gi. 



A dependence of the structure function gi{x,Q^) given by Eqs ([I|) - (|p on the 
parameter C in Eqs (^) and (^ is illustrated in Fig. 2 where we plot the asymmetries 
Ai{x) for the proton at the measured and for < 1 GeV^. We expect the VMD 
contribution to be dominant there. This cut selected the SMC data [0, ^ at low values 
of X and the SLAC E143 measurements at 16.2 GeV incident electron energy at 
higher a;.0 To obtain predictions for the asymmetry Ai rather than for gi we used the 
model |Tl| for the F2 and two different parametrisations [0, ^ for R, as employed in 
the E143 and SMC analyses respectively. 



The statistical accuracy of the SMC data is too poor to constraint the value of 
the coefficient C, i.e. of the VMD-type nonperturbative contribution to the struc- 
ture function gi{x,Q'^) at low values of Q^. The SLAC E143 data apparently prefer 
a small negative value of C. The model prediction without VMD contribution (C=0) 
is systematically higher than the EMS measurements. The fact that the data prefer 
negative value of the VMD contribution is consistent with the results obtained from 



the phenomenological analysis of the sum-rules [ 25 



Similar analysis performed for the neutron and deuteron structure functions, g^ 
and gf, where in the former case data cover narrower kinematic interval and in the 
latter the statistics at low x is substantially poorer, turned out to be inconclusive. 



5. Summary and conclusions 

We have analysed the recent gi{x, Q"^) measurements at low values of x and within 

a formalism based on unintegrated spin dependent parton distributions incorporating 

the leading order Altarelli-Parisi evolution and the double ln^(l/x) resummation at 

low X. A VMD-type nonperturbative part was also included since low values of x in 

"'The E143 measured Ai for < 1 GcV^ also at 9.7 GeV incident electron energy. For these data 
4< ^ 10 GeV^, i.e. above the resonance region but too small for our model to be applicable. 
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the measurements correlate with low values of Q^. The ln^(l/x) effects are not yet 
important in the kinematic range of the fixed target experiments but the formalism 
based on unintegrated parton distributions, summarised by Eq.(|^), is very suitable for 
extrapolating gi to the region of low . The model reproduces a general trend in the 
data for the proton. The statistical accuracy of the SMC measurements taken at lowest 
values X, x >0. 00006, and of the Q^, > 0.01 GeV^, is however too poor to constraint 
the VMD contribution. A more accurate data from the SLAG E143 experiment, where 
X > 0.02 and > 0.5 GeV^ seem to prefer a nonzero and negative contribution of the 
VMD to Qi of the proton. 
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Figure 1: Summary of the xgi measurements for the proton and for the deuteron as 
a function of x at the measured obtained with different experiments. The inserted 
figures show the SMC data for which < 1 GeV^. Errors are statistical The 
curves, calculated at x and values at the SMC measurements result from the model 
described in Sec. 3. 
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Figure 2: The asymmetry Ai for the proton as a function of x at the measured 
(marked above the x axis), obtained by the SMC |, i and SLAG E143 g (at 16.2 
GeV incident energy). Errors are statistical. Curves are calculated according to Eqs 
(HI) ^ (H) assuming different values of C in Eqs (|^) and (H). 
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